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ABSTRACT 


The sea state bias (SSB) has become the dominant 
source of error in satellite altimetry. The operational SSB 
correction models are two-dimensional (2-D) nonparametric 
models based on the wind speed (U) and the significant 
wave height (SWH) that can be directly measured by the 
altimeters. This paper estimates the sea state bias of HY-2A 
radar altimeter using a three-dimensional (3-D) 
nonparametric model based on SWH from HY_2A interim 
geophysical dataset records (IGDR), U and the mean wave 
period (MWP) from the European Centre for Medium- 
Range Weather Forecasts (ECMWF) reanalysis project 
ERA-Interim. The 3-D SSB estimates can increase the 
explained variance by 1.72 cm’, or 1.31 cm RMS relative to 
the traditional 2-D SSB estimates based on U and SWH. 


Index Terms— HY-2A, radar altimeter, sea state bias, 
crossover differences, 3-D nonparametric model 


1. INTRODUCTION 


Sea level rise mainly caused by global warming has attracted 
more and more people’s attention. One of the main purposes 
of radar altimeters is to measure the mean sea level. 
However, the sea level measured by radar altimeters is lower 
than the true sea level because wave troughs are better 
reflectors of radar energy than wave crests. This effect is 
called the electromagnetic bias (EM bias). There are two 
other sea-state-related biases: the skewness and tracker 
biases[1]. The EM bias together with the skewness bias and 
the tracker bias is generally called the sea state bias (SSB). 
Thanks to the improvement of the precise orbit 
determination technology and other geophysical corrections, 
the sea state bias has become the dominant source of error in 
satellite altimetry. 

The operational SSB correction models are two- 
dimensional (2-D) empirical models based on the wind 


speed (U) and the significant wave height (SWH) that can be 
directly measured by the altimeters. However, these 2-D 
SSB models can not entirely parameterize the range bias 
variability and the SSB uncertainty may be lowered through 
improved SSB models by including additional measurable or 
predictable correlatives [2]. This paper estimates the sea 
state bias of HY-2A radar altimeter by using a 3-D 
nonparametric model based on U, SWH and MWP. The 
model is developed using the crossover sea surface height 
(SSH) differences. Compared with the SSB corrections in 
HY-2A IGDR and the traditional 2-D SSB estimates, the 3- 
D SSB estimates can greatly reduce the standard deviation 
(STD) of the crossover SSH differences and the STD of the 
sea level anomalies (SLAs). The enhancement from 2-D to 
3-D SSB estimate is of great significance for improving the 
precision of HY-2A altimeter product. 


2. METHODOLOGY 


In this paper, we developed the 3-D SSB model using 
crossover SSH differences from HY-2A altimeter. 

The sea surface height measurement uncorrected for the 
SSB is denoted as SSH . SSH includes the geoid (h, ), the 
dynamic topography (y), the SSB and some measurement 
noise (@): 

SSH =h,+7+SSB+@ (1) 

At crossover points, the geoid signal can be eliminated by 
forming differences between two SSH measurements taken 


on the ascending arc (index 1) and the descending arc (index 
2) respectively. 


SSH, — SSH, = (SSB, —SSB,)+(7,-7,)+(@,-@,) (2) 


Assume (77, —77,)+(@,—@,) can be considered as a single 
zero-mean noise termé. The observation equation (2) for 
the SSH differences can be written as: 


SSH, - SSH, = SSB, - SSB, + £ (3) 


Assume the SSB is a function g ofx =(U, SWH, MWP): 
SSB = 9(x) (4) 
Noting y = SSH, - SSH, , the observation equation (3) can be 
written as: 
Y =GA(X,)- P(X) +é (5) 
Equation (5) can be solved by using a nonparametric 


estimation method based on kernel smoothing introduced by 
Gaspar [3, 4]. 
3. DATA 
The input parameters of the 3-D SSB model developed 
are U, SWH and MWP. U and SWH usually are obtained 


from the altimeter observations. However, the wind speed in 
HY-2A IGDR is not good, so we used the wind speed from 
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ECMWF reanalysis project ERA-Interim. Fig.1 gives U and 
SWH from ECMWF and Jason-2 GDR respectively. Jason- 
2 altimeter has been fully calibrated. Both ECMWF derived 
U and SWH have high accuracy compared to Jason-2 GDR 
data. Fig.2 gives U and SWH from ECMWF and HY-2A 
IGDR respectively. The root mean square error (RMSE) of 
the differences between ECMWF derived SWH and HY-2A 
IGDR derived SWH is only 0.365m. Therefore, the SWH 
from HY-2A IGDR can be used to develop the SSB model. 
The RMSE of the differences between ECMWF derived U 
and HY-2A IGDR derived U is larger than 5m/s (The 
backscattering coefficient in HY-2A IGDR has not been 
fully calibrated, after the backscattering coefficient is fully 
calibrated, the derived U will be better), so we use ECMWF 
derived U to estimate the SSB of HY-2A altimeter. 
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Fig.1 U and SWH comparison between ECMWF and Jason-2 altimeter (left: U; right: SWH) 
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Fig.2 U and SWH comparison between ECMWF and HY-2A altimeter (left: U; right: SWH) 


3.1 ERA-Interim reanalysis 


ERA-Interim is the third global reanalysis project 
produced by the European Centre for Medium-Range 
Weather Forecasts (ECMWF). ERA-Interim provides 
reanalysis data from 1979 to the present time and updates in 
near real time. The objective of the ERA-Interim project is 
to prepare for a new global reanalysis project to replace 
ERA-40 gradually. Compared with ERA-40, ERA-Interim 
has some improvements in data assimilation, forecast model 
and observing system. A variety of data in uniform 


latitude/longitude grids can be obtained from ERA-Interim 
datasets. The mean wave period and the wind speed we used 
are projected on the grid of 0. 25° x 0.25°. 


3.2 HY-2A altimeter 


HY-2A satellite was launched on 16 August, 2011. HY- 
2A is a polar-orbit satellite with an orbit height of 971 km. 
The instruments on HY-2A include a dual-band (13.58 and 
5.25 GHz) radar altimeter (RA), a Ku-band radar 
scatterometer (SCAT) and a microwave imager (MWD). The 
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RA is equipped with a three-channel (18.7, 23.8 and 37 GHz) 
microwave radiometer for wet troposphere delay correction. 

In this paper, we used the cycle 42-56 of the HY-2A 1-Hz 
IGDR data to develop the SSB models. 


4. RESULTS 


4.1 HY-2A SSB estimates 


The 2-D SSB estimates are shown in Fig.3 in the form of 
2-D grids with contours given in cm. As is shown in Fig.3, 
the magnitude of the SSB generally increases with the 
increase of SWH. For a given SWH, the magnitude of the 
SSB first increases with U and then decreases at higher wind 
speeds. 
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Fig.3 HY-2A SSB estimates (in cm) obtained using a 2-D SSB 
model based on U and SWH. Shaded areas represent data density: 
(dark gray) when it is less than 5 samples, (red) when it is between 
5 between 30, (yellow) when it is between 30 and 100, and (white) 
when it is larger than 100 samples. 


The 3-D SSB model based on U, SWH and MWP is 
presented as a 3-D lookup table and the ranges of U, SWH 
and MWP are respectively 0-30m/s, 0-12m and 0-18s. In 
order to better represent the 3-D SSB model, we present the 
correction values in the form of 2-D arrays by fixing a third 
parameter. Contours of the SSB values are displayed in cm 
and they are shown in Fig. 4. The fixed values are near the 
mean values to ensure that there are enough measurements. 
In our study, the three fixed value we choose are MWP=9s, 
U=7.5m and SWH=2.5m respectively. Fig.4 (a) agrees 
closely with Fig.3. When U and MWP are fixed, the 
magnitude of the SSB is an increasing function of SWH. In 
contrast, as can be seen in Fig.4 (b), the magnitude of the 
SSB is a decreasing function of MWP when U and SWH are 
fixed. In the data-rich region, there is an approximately 4 cm 
variation versus MWP, with the magnitude of the SSB 
decreasing as the wave period increases. As can be observed 
in Fig.4 (c), when SWH is fixed, variations with MWP are 
more obvious than that with U in data-rich region. For a 
given U, the magnitude of the SSB is a decreasing function 
of MWP. 


4.2 Assessment of the 3-D SSB estimates 


We used the explained variance to evaluate the 3-D SSB 
model. The explained variance is the reduction of the 
variance of the crossover SSH differences obtained after 
applying the given SSB correction and it is a main criterion 
to assess the quality of a SSB model. We first used the 
global data to calculate the explained variances. The 
variances explained by 3-D SSB estimates, 2-D SSB 
estimates and the SSB corrections in HY-2A IGDR are 
29.95, 28.23 and 26.47 cm”. 
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Fig.4 HY-2A SSB estimates (in cm) obtained using a 3-D SSB 
model based on U, SWH and MWP. This is shown using three 2-D 
arrays with the respective third variable held constant. The fixed 
values are (a) MWP=9s, (b) U=7.5m/s, (c) SWH=2.5m 
respectively. Shaded areas represent data density: (dark gray) 
holding no data, (red) when it is between 1 between 5, (yellow) 
when it is between 5 and 10, and (white) when it is larger than 10 
samples. 


We then binned the crossover data into 10° latitude bands 
and computed the explained variance in each band. The 
explained variance in comparison to the 1D (-3.8%SWH) 
SSB model is given in Fig.5. 3-D SSB estimates are better at 
most latitude regions. 
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Fig.5 Latitude distribution of the variance explained by 3D SSB 

estimate and 2D SSB estimate respectively in comparison to a 1D 

(-3.8%SWH) benchmark. The variance explained by the SSB 


corrections in HY-2A IGDR data is also given. 


Crossover SSH differences are the main tool to analyze 
the whole altimetry system performance. They allow us to 
analyze the SSH consistency between ascending and 
descending passes [5]. Fig.6 gives the standard deviation 
(STD) of the crossover SSH differences after using different 
SSB corrections. Results are obtained cycle per cycle. 
Obviously, the 3-D SSB estimates can greatly reduce the 
STD of the crossover SSH differences. 
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Fig.6 The standard deviation (STD) of the crossover SSH 
differences after using 3-D SSB estimates, 2-D SSB estimates and 
the SSB corrections in HY-2A IGDR data respectively. 


Along track analyze is also used to assess the altimeter 
system performance by computing the sea level anomalies 
(SLAs). SLAs are computed along track by subtracting the 
mean sea surface (MSS) from the instantaneous SSH 
measurements. The SLAs variance gives an estimate of the 
errors of the system [6]. Fig.7 gives the STD of SLAs after 
using different SSB corrections. The STD of SLAs can be 


reduced after replacing the traditional 2-D SSB estimates 
with the 3-D SSB estimates. 
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Fig.7 The standard deviation (STD) of the along-track sea level 
anomalies (SLAs) after using 3-D SSB estimates, 2-D SSB 
estimates and the SSB corrections in HY-2A IGDR data 
respectively. 


5. CONCLUSIONS 


This paper estimates the SSB of HY-2A radar altimeter 
using a 3-D nonparametric model. Crossover SSH 
differences are used to develop the SSB model. Compared 
with the traditional 2-D SSB estimates, the 3-D SSB 
estimates can greatly reduce the STD of the crossover SSH 
differences and SLAs. Therefore, the 3-D SSB estimates can 
improve the precision of HY-2A altimeter product. 
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